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We investigate the holographic dark energy scenario with a varying gravitational constant, in 
flat and non-flat background geometry. We extract the exact differential equations determining 
the evolution of the dark energy density-parameter, which include G-variation correction terms. 
Performing a low-redshift expansion of the dark energy equation of state, we provide the involved 
parameters as functions of the current density parameters, of the holographic dark energy constant 
and of the G-variation. 
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I. INTRODUCTION 



vacuum energy as holographic dark energy: 
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Recent cosmological observations obtained by SNe la 
[H, WMAP 0, SDSS 3 and X-ray [| indicate that 
the universe experiences an accelerated expansion. Al- 
though the simplest way to explain this behavior is the 
consideration of a cosmological constant [sj] , the two rel- 
evant problems (namely the "fine-tuning" and the "co- 
incidence" one) led to the dark energy paradigm. The 
dynamical nature of dark energy, at least in an effective 
level, can originate from various fields, such is a canon- 
ical scalar field (quintessence) [6|, a phantom field, that 
is a scalar field with a negative sign of the kinetic term 
Q, or the combination of quintessence and phantom in 
a unified model named quintom [Sj]. 

Although going beyond the above effective description 
requires a deeper understanding of the underlying the- 
ory of quantum gravity Q unknown at present, physi- 
cists can still make some attempts to probe the nature 
of dark energy according to some basic quantum gravita- 
tional principles. An example of such a paradigm is the 
holographic dark energy scenario, constructed in the light 
of the holographic principle [13, El, [13, (although the 
recent developments in Horawa gravity could offer a dark 
energy candidate with perhaps better quantum gravita- 
tional foundations 14]). Its framework is the black hole 
thermodynamics and the connection (known from 
AdS/CFT correspondence) of the UV cut-of of a quan- 
tum field theory, which gives rise to the vacuum energy, 
with the largest distance of the theory Thus, de- 

termining an appropriate quantity L to serve as an IR 
cut-off, imposing the constraint that the total vacuum 
energy in the corresponding maximum volume must not 
be greater than the mass of a black hole of the same size, 
and saturating the inequality, one identifies the acquired 
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with G the Newton's gravitational constant and c a con- 
stant. The holographic dark energy scenario has been 
tested and constrained by various astronomical observa- 
tions [M 113, H [iljpojl and it has been extended to 
various frameworks 2Vi 2^ . [23j . 



Until now, in all the investigated holographic dark en- 
ergy models a constant Newton's "constant" G has been 
considered. However, there are significant indications 
that G can by varying, be ing a function of time or equiv- 
alcntly of the scale factor [2J]. In par ticular, observations 
of Hulse- Taylor binary pulsar [25ll26j] , helio-seismological 
data 27], Type la supernova observations [l| and aster- 
eoseismological data from the pulsating white dwarf star 



G117-B15A \m lead to 



G/G 



< 



4.10 X IQ-^^yr-^, 



for 



z < 3.5 [30|. Additionally, a varying G has some theoret- 
ical advantages too, alleviating the dark matter problem 
[31I ] , the cosmic coincidence problem (32I] and the discrep- 
ancies in Hubble parameter value [33|. 

There have been many proposals in the literature at- 
tempting to theoretically justified a varying gravitational 
constant, despite the lack of a full, underlying quantum 
gravity theory. Starting with the simple but pioneering 
work of Dirac [3^ , the varying behavior in Kaluza-Klein 
theory was associated with a scalar field appearing in the 
metric component corresponding to the 5-th dimension 
[35! and its size variation ^3^. An alternative approach 
arises from Brans-Dickc framework (37j , where the gravi- 
tational constant is replaced by a scalar field coupling to 
gravity through a new parameter, and it has been gen- 
eralized to various forms of scalar- tensor theories [38| . 
leading to a considerably broader range of variable-G 
theories. In addition, justification of a varying Newton's 
constant has been established with the use of conformal 
invariance and its induced local transformations [3^ . Fi- 
nally, a varying G can arise perturbatively through a 
semiclassical treatment of Hilbert- Einstein action [40| . 
non-perturbatively through quantum-gravitational ap- 



2 



proaches within the "Hilbert-Einstein truncation" [41l |. 
or through gravitational holography \^ . 

In this work we are interested in investigating the 
holographic dark energy paradigm allowing for a vary- 
ing gravitational constant, and extracting the corre- 
sponding corrections to the dark energy equation-of- 
state parameter. In order to remain general and ex- 
plore the pure varying-G effects in a model-independent 
way, we do not use explicitly any additional, geometri- 
cal or quintessence-like, scalar field, considering just the 
Hilbert-Einstein action in an affective level, as it arises 
from gravitational holography [i^, |4^. In other words, 
we effectively focus on the dark energy and dark mat- 
ter sectors without examining explicitly the mechanism 
of G-variation, which value is considered as an input 
fixed by observations. Additionally, generality requires 
to perform our study in fiat and non-flat FRW universe. 
The plan of the work is as follows: In section |TT] we con- 
struct the holographic dark energy scenario with a vary- 
ing Newton's constant and we extract the differential 
equations that determine the evolution of dark energy 
density-parameter. In section IIIII we use these expres- 
sions in order to calculate the corrections to the dark 
energy equation-of-state parameter at low redshifts. Fi- 
nally, in section HV] we summarize our results. 



II. HOLOGRAPHIC DARK ENERGY WITH 
VARYING GRAVITATIONAL CONSTANT 

A. Flat FRW geometry 

Let us construct holographic dark energy scenario al- 
lowing for a varying Newton's constant G. The space- 
time geometry will be a fiat Robertson- Walker: 



with a{t) the scale factor and t the comoving time. As 
usual, the first Friedmann equation reads: 



8ttG 



(Pm + pa) , 



with H the Hubble parameter, p„j — where pm 

and PA stand respectively for matter and dark energy 
densities and the index marks the present value of a 
quantity. Furthermore, we will use the density parameter 
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PA, which, imposing explicitly the holographic 



In the following we will use In a as an independent 
variable. Thus, denoting by dot the time-derivative and 
by prime the derivative with respect to In a, for every 
quantity F we acquire F — F'H . Differentiating @ 
using ([5]), and observing that Rh = HRh — 1, we obtain: 



2 r 



- 1 - 



H 



(6) 



Until now, the varying behavior of G has not become 
manifested. However, the next step is to eliminate H . 
This can be obtained by differentiating Friedman equa- 
tion ([3]), leading to 



H 

2— = -3 + n, 




^a) , (7) 



where G is considered to be a function of In a. In the 
extraction of this relation we have additionally used the 
auxiliary expression 



(8) 



which arises from differentiation of ([T]). Therefore, sub- 
stituting ([7]) into ([6]) we finally obtain: 




G ■ 



(9) 



The first term is the usual holographic dark energy dif- 
ferential equation [T3| . The second term is the correction 
arising from the varying nature of G. Note that G'/G is 
a pure number as expected. 

Finally, for completeness, we present the general solu- 
tion for arbitrary c and G'/G = Ac which in an implicit 
form reads 



(2) 1^ 




ln(l - 



2 + c(l- Ag) 
81n[c(l-- AG) + 2x/ni] 



c(Ag-1) c(Ag-1)[c2(Ag-1)2-4] 



(10) 



(3) The constant xq can be straightforwardly calculated if we 
determine aq and fij^ today (for example choosing — 1 



nature of dark energy according to relation ([1]) , becomes 



Xo is equal to the left hand side with Vt^ replaced by ^2^)- 
Clearly, for Ag = and c—1, expression (fTU]) coincides 
with that of UM- 



B. Non-flat FRW geometry 



Oa = 



(4) 



Finally, in the case of a flat universe, the best choice for 
the definition of L is to identify it with the future event 
horizon [H, HI, S, ill , that is L = Rh{a) with 



Rh{a) 



dt' 



da' 



(5) 



In this subsection we generalize the aforementioned 
analysis in the case of a general FRW universe with line 
element 



ds^ = -df + a^{t) 



dr^ 



1 ~ kr^ 



(11) 



in comoving coordinates {t,r,9,ip), where k denotes the 
spacial curvature with k = —1,0,1 corresponding to 
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open, flat and closed universe respectively. In this case, 
the first Fricdmann equation writes: 



2 k SnG ( 
a-' 



PA 



(12) 



According to the formulation of holographic dark en- 
ergy in non-flat geometry, the cosmological length L in 
(HI) is considered to be [2l| : 



a{t) 



smn 



klRh 



(13) 



where 



1 



sin 2/ fc = +1, 

smn{y/\k\y) = fc = 0, 

sinhy fc = — 1. 

A straightforward calculation leads to 



(14) 



(15) 



where 



cosy fc — +1, 
cosn(-\/|fc|y) = 1 fc = 0, 

coshy fc = — 1. 



(16) 



Repeating the procedure of the previous sub-section 
and differentiating (g]) using and ([TS]) we obtain: 

On the other hand, differentiating Fricdmann equation 
(fT2ll we finally obtain 



= — 3 — fife + r^A 



3/2 



■ cosn 



\k\Rh 



+ {1 + ilk -Ha) 



G" 



(18) 



where we have introduced the curvature density param- 
eter rife = ^^ J^yi ■ Therefore, substituting into (fTT]) 
we result to 



l + D,k—i^A + —^—— cosn 



|fc|i?h 



G' 



-nAii + iik-nA)—. (19) 

Expression ^9]) provides the correction to holographic 
dark energy differential equation in non-flat universe, due 
to the varying nature of G. Clearly, for fc = (and thus 
= 0) it leads to 



III. COSMOLOGICAL IMPLICATIONS 

Since we have extracted the expressions for Q'^, we 
can calculate w{z) for small redshifts z, performing the 
standard expansions of the literature. In particular, since 
PA ^ a~^^^+"') we acquire 



In PA = In Pa 



, d\npA 1 d^lnpA 



dlna^""" ' 2d(lna)2 



(lna)2 + ..., (20) 



where the derivatives are taken at the present time oo = 1 
(and thus at J^a = ^a)- Then, w(lna) is given as 



w{\na) = -1 - - 



tiln PA 1 cPjnpA , 
rflna 2d(lna)2 



(21) 



up to second order. SincepA = 3H^Q\/{8ttG) = 
^Apm/^m = Pmof^A/(l + ^^fc — ^A)a^^ , the derivatives 
are easily computed using the obtained expressions for 
^I'a- In addition, we can straightforwardly calculate w{z), 
replacing In a = — ln(l-|-2;) ~ — z, which is valid for small 
redshifts, defining 



w{z) 



d In PA ' 
din a 





(P In PA 




_d(lna)2_ 



Z = Wq + WlZ. 

(22) 

The role of G- variation will be expressed through the 
pure number G' /G = Aq, which will be extracted from 
observations. In particular, observations of Hulse- Taylor 
binary pulsar B1913 -|- 16 lead to the estimation G/G ^ 



2±4 X 10 yr ^ [25|, |26|, while helio-seismological data 
provide the bound -1.6 x lO^^^^r-i < G/G < 
Similarly, Type la supernova observations [l'] give the 
best upper bound of the variation of G as — 10~^^yr~^ < 
^ < at redshifts z : 



0.5 [28|, while astereoseismological 
data from the pulsating white dwarf star G117-B15A lead 

to g < 4.10 X 10-"yr-i See also [sS] for various 

bounds on G/G from observational data, noting that all 
these measurements are valid at relatively low redshifts, 
i.e z < 3.5. 

Since the limits in G- variation are given for G/G in 
units yr~^, and since G/G = HG' /G, we can estimate 
Ag substituting the value of H in yr^^ . In the fol- 
lowing we will use g < 4.10 x lO^^^j/r^^. Thus, in- 
serting an average estimation for the Hubble parame- 
ter « (i?) « 6 X 10-"?/r-i [Hi, we obtain that 
< |Ag| < 0.07. Clearly, this estimation is vahd at 
low redshifts, since only in this range the measurements 
of G/G and the estimation of the average (H) are valid. 
However, the restriction to this range is consistent with 
the z-expansion of w considered above. 
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A. Flat FRW geometry 

In this case ft'j^ is given by ([9]), and the aforementioned 
procedure leads to 



Wo 



Wl 




(23) 



(24) 



These expressions provide wq and wi , for the hofographic 
dark energy with varying G, in a flat universe. Obviously, 
when Ag — 0, they coincide with those of [isj . 

In general, apart from the relevant uncertainty in 17^ 
measurements, we face the problem of the uncertainty 
in the constant c. In particular, observational data from 
type la supernovae give the best-fit value c = 0.21 within 
1-cr error range jl^ , while those from the X-ray gas mass 
fraction of galaxy clusters lead to c = 0.61 within 1-a 
[l3|. Similarly, combining data from type la supernovae. 
Cosmic Microwave Background radiation and large scale 
structure give the best-fit value c = 0.91 within l-cr p^ . 



while combining data from type la supernovae. X-ray 
gas and Baryon Acoustic Oscillation lead to c = 0.73 
as a best- fit value within 1-cr p^ . However, expressions 
(ESI I (EH) provide the pure change due to the variation of 
gravitational constant for given c and For example, 
and in order to compare with the corresponding result 
of [13], imposing $1° « 0.73 and c = 1, and using < 
|Ag| < 0.07 we obtain: 



Wo 



-0.903: 



Wl = 0.1041 



-0.023 
-0.023 
4-0.0025 
-0.0025' 



(25) 



where we have neglected uncertainties other than G- 
variation. Finally, note that the wp-variation due to Aq 
is absolute, that is it does not depend on c and fi^ , while 
that of Wl does depend on these parameters. However, 
the relative variations of wo,wi do depend on the c- value, 
and they are smaller for smaller c. 



B. Non-flat FRW geometry 

In this case il^ is given by and the aforemen- 

tioned procedure leads to 



Wo 



Wl 



k\Rho 



6c 



1 + nl 



2 



Ag 
3 

cosn 



k\Rha 



(1 



(26) 




(27) 



In these expressions, ftl is the present day value of the 
curvature density parameter, and we have defined 

{sin^y k = +1, 

^ fc-0, (28) 
— sinh" y k — —1. 

Finally, i?^o ^-nd oq are the present values of the corre- 
sponding quantities. Clearly, for A: = 0, that is for a flat 
geometry, (PB)) . (|77|) coincide with ([^ . ([M)) respectively. 

As we observe, expressions (|26t . (|27)) . apart from the 
present values of the parameters ri°, fl^ contain ao 
and the value of Rho at present. This last term is 
present in a non-flat universe, and it is a "non-local" 
quantity which has to be calculated by an integration 
(see relations (IT51) and (O). However, making use of 



the holographic nature of dark energy, we can over- 
come this difficulty. Indeed, from (jH) we obtain that 
Lq — c/ (Hoy^Q'j^), with Hq the present value of the Hub- 
ble parameter. On the other hand, from (|13p we acquire 

Rho/ttQ ~ } sinn~"'"(-\/[fc[Lo/ao). Therefore, we con- 
v 1^ 

elude that 

Rm ^ ■ -if \ 

= — ^smn — r= — 

ao ^\ {aoHoy^J 

^^smn-MV (29) 

a relation which proves very useful. Substituting into 
([26|) ■ ([27]) we finally obtain the simple expressions: 
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(30) 



Wi 



6c 



'm (i + nl- nl) Ag. (31 



Note that wo,wi depend eventually only on O^, ftf,, c and 
of course Aq. Similarly to the previous subsection, in 
order to give a representative estimation and neglecting 
uncertainties of other quantities apart from G-variation, 
we use c = 1, « 0.73, « 0.02, < |Ag| < 0.07, 
obtaining: 

wo = -0.895t°°23 

«;i= 0.111™. (32) 

Finally, we mention that the relative variations of wo,wi 
depend on the c-value, and they are smaller for smaller 
c. 

IV. CONCLUSIONS 

In this work we have investigated the holographic dark 
energy scenario with a varying gravitational constant, go- 
ing beyond the simple scenarios of [4^. Imposing flat 
and non-flat background geometry we have extracted 
the exact differential equations that determine the evo- 
lution of the dark energy density-parameter, where the 
G-variation appears as a coefhcient in additional terms. 
Thus, performing a low-redshift expansion of the dark 
energy equation-of-state parameter w{z) wq + wiz, we 
provide wo,wi as functions of fi^, fi^, of the holographic 



dark energy constant c, and of the G-variation Ag (ex- 
pressions (|30p . (PT|) '). As expected, the variation of the 
gravitational constant increases the variation of w{z). 

In the aforementioned analysis, the G-variation has 
been considered as a constant quantity at the cosmo- 
logical epoch of interest, that is at low redshifts, as it 
is measured in observations with satisfactory accuracy 
25, 26, 27, 28, 29, 30]. A step forward would be to 
consider possible G(2;)-parametrizations 0, \^ and ex- 
tract their effect on w{z). However, such parametriza- 
tions have a significant amount of arbitrariness, since the 
present observational data do not allow for such a reso- 
lution, and thus we have not performed this extension in 
the present work. 

Finally, we mention that in general, the possible un- 
certainty of the constant c can have a larger effect on 
w{z) than that of G-variation. In the above investiga- 
tion we have just provided the complete expressions, in- 
cluding the correction terms due to the variation of the 
gravitational constant. One could proceed to a combined 
observational constraint analysis, allowing for variations 
and uncertainties in all parameters, as it was part ially 
performed in the specific Brans-Dicke framework in [47| . 
This extended examination, with not-guaranteed results 
due to complexity, is under current investigation and it 
is left for a future publication. 
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